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Abstract 

We discuss the parameter space of the Inert Doublet Model, a two Higgs doublet model with a dark 
matter candidate. An extensive set of theoretical and experimental constraints on this model is considered, 
where both collider as well as astroparticle data limits, the latter in the form of dark matter relic density as 
well as direct detection, are taken into account. We discuss the effects of these constraints on the parameter 
space of the model. In particular, we do not require the IDM to provide the full dark matter content of 
the universe, which opens up additional regions in the parameter space accessible at collider experiments. 
The combination of all constraints leads to a relatively strong mass degeneracy in the dark scalar sector 
for masses > 200 GeV, and to a minimal scale ~ 45 GeV for the dark scalar masses. We also observe a 
stringent mass hierarchy > M^. We propose benchmark points and benchmark planes for dark scalar 
pair-production for the current LHC run being in compliance with all theoretical as well as experimental 
bounds. 
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1 Introduction 


1. INTRODUCTION 

The Inert Doublet model (IDM) is one of the most straightforward extensions of the Standard 
Model (SM) [IHS]. It is a special type of Two Higgs Doublet Model (2HDM). In the scalar sec¬ 
tor it contains, in addition to the SU(2) doublet ips responsible for the spontaneous breaking of 
electroweak symmetry, a second scalar doublet (I)d with vanishing vacuum expectation value {vev). 
This second doublet is not involved in mass generation via spontaneous symmetry breaking and 
does not couple to fermions (in this respect the model corresponds to a Type I 2HDM.). In this 
model a Z 2 symmetry, which we label D-symmetry below, with the discrete transformations defined 
as : 


(pD ~(pDi (ps P’S-! SM —)• SM, (1) 

is respected by the Lagrangian and the vacuum. 

The ps doublet plays the same role as the corresponding doublet in the SM, providing the 
SM-like Higgs particle. This doublet is even under the D symmetry, while the second doublet, pD, 
is D-odd. This so called inert or dark doublet contains 4 scalars, two charged and two neutral ones, 
with the lightest neutral scalar being a natural DM candidate. This model was studied in order 
to provide a heavy Higgs boson [3] as well as a lighter Higgs boson, to be produced at LHC [2] . It 
was considered as a model with a “perfect example” of a weakly interacting massive particle [M]. 
It leads to an interesting pattern of the Universe evolution, towards the Inert phase as given by 
the IDM, with one, two or three phase transitions [8j. Furthermore, the IDM can provide a strong 
first-order phase transition [MS], which is one of the Sakharov conditions needed to generate a 
baryon asymmetry of the Universe. After the discovery of a SM-like Higgs particle in 2012, many 
analyses have been performed for the IDM, using Higgs collider data, as well as astrophysical 
measurements, see e.g. In addition, proposals were made how to search for dark scalars 

at the LHC in the di- or multilepton channel [IS1E21E3]. 

Recently, also the important issue of vacuum (meta)stability in the IDM has been discussed, 
and it was found that additional (possibly heavy) scalars can have a strong impact on it [23II27] 
Especially facing the LHC run H, the determination of the regions of parameter space which survive 
all current theoretical and experimental constraints is indispensible. Therefore, we here provide a 
survey of the IDM parameter space, after all constraints are taken into account. The main sources 
of these stem from run I LHC Higgs data as well as dark matter relic density (Planck) and direct 
dark matter search at LUX. A hrst preliminary preview on some of our results was presented in a 
conference proceeding m- 


^ Recent analyses for models which extend the IDM by an additional singlet have been performed in |18H21 |. 

^ Similar solntions can be fonnd in a simple singlet extension of the SM Higgs sector, cf. e.g. |28[l29j and references 
therein. 
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2 The model 


2. THE MODEL 


The D-symmetric potential of the IDM has the following form: 


( 2 ) 


with all parameters real (see e.g. [ 8 ]). In a 2HDM with this potential different vacua can exist, 
e.g. a mixed one with (^ 5 ) 7 ^ 0 , (<?i>D) 7 ^ 0 , an inertlike vacuum with {(ps) = 0 , (<^d) / 0 , or even 
a charge breaking vacuum (see [3 0 EH EU). The decomposition around the vacuum state in the 
IDM is given by 

where v = 246 GeV denotes the vacuum expectation value. 

The first doublet, contains the SM-like Higgs boson h, with the mass 

Ml = Aiu^ = (4) 


which, according to the LHC data, is close to 125 GeV [33] . 

The second doublet, (pD, consists of four dark (inert) scalars H, A, H^, which do not couple to 
fermions at the tree-level. Their masses are given as follows: 

\ - ^22) , 

M\ = M|^± -h ^ (A 4 - As) = ^(A345?^^ - ^ 22 ), 

Mfj = Mfj± + — (A 4 -|- As) = -(A 34 sU^ — ?Ti 22 )- (5) 

Due to an exact D symmetry the lightest neutral scalar H (or A) is stable and thereby may serve 
as a good dark matter candidate.^ We take H to be the DM candidate and so Mh < Ma, 

(this choice implies As < 0 , A 4 S = A 4 -|- As < 0 ). Note that, as the second doublet does not couple 

to fermions, one cannot ascribe a CP-property to the scalars of the dark sector. Therefore, in 
contrast to general two Higgs doublet models, in our scenario the two neutral scalars H and A 
can be treated on equal footing, and in fact the corresponding couplings correspond to each other 
using the replacement As o — As, cf. Appendix 
The parameters 


A345 — A3 -|- A4 -|- As, A34S — A3 -|- A4 — As 


( 6 ) 


® Charged DM has been strongly limited by astrophysical analyses [34] • 
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3 Constraints 


3.1 Theoretical constraints 


are related to the triple and quartic coupling between the SM-like Higgs h and the DM candidate 
H or the scalar A, respectively, while A 3 describes the h interaction with the charged scalars 
The parameter A 2 describes the quartic self-couplings of dark particles. Note that an equivalent 
choice would be to take A as the lightest neutral dark scalar and therefore as the dark matter 
candidate, which implies A 5 > 0. In this case, the discussion below remains valid by replacing 

As ^^ —As, A 345 <—A 34 S. 

The parameter m |2 sets here a common mass scale for the dark particles. For a list of all relevant 
Feynman rules for this model, see Appendix 

After electroweak symmetry breaking, the model has seven free parameters. Since the val¬ 
ues of V and the SM-like Higgs mass are fixed by experimental data, there are in total 5 free 
parameters defining the scalar sector of the IDM. Typical choices are either physical parameters 
{Mh, Ma, Mjj±, A 2 , A 345 ) or potential parameters (r?T, 22 , A 2 , A 3 , A 4 , A 5 ). We here choose to work in 
the physical basis. 


3. CONSTRAINTS 

As has been widely discussed in the literature, the IDM is subject to numerous constraints, 
both from theoretical conditions as well as from experimental results. We here briefly remind the 
reader of these constraints and refer to the literature [a nn [la Ea I35H15] for further details. 


3.1. Theoretical constraints 

• First, we require the potential to be bounded from below, such that there is no field configu¬ 
ration for which V —)• —00. This directly leads to the following conditions on the couplings 

m 

Ai > 0, A 2 > 0, As -b -\/A 1 A 2 > 0, A 345 -b \/A 1 A 2 > 0, (7) 

with A 345 as defined above. The above relations hold on tree level and neglect higher or¬ 
der contributions which could (de)stabilize the electroweak vacuum. For a more detailed 
discussion of such effects, cf. e.g. [25l [26] . 

• We also require the scalar 2 —)• 2 scattering matrix^ to be unitary. We do this by the means 
of perturbative unitarity [48| . where we consider all scattering matrices for scalars with 
specific hypercharge and isospin and require the eigenvalues to be |Lj| < Svr. This follows 
the generic prescription as e.g. given in [l9], where we make use of the routine implemented 
in the Two Higgs Doublet Model Calculator (2HDMC) tool [50]. 


We consider electroweak gauge boson scattering in the high energy limit, where the corresponding Goldstone modes 
are used in accordance with the equivalence theorem m- 
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3 Constraints 


3.2 Experimental constraints 


• Furthermore, we require all quartic Higgs coupling to be perturbative, i.e. to take absolute 
values < An. We apply this limit both to the coupling parameters in the potential, Aj, as 
well as the couplings stemming from vertex Feynman rules. ^ 

• In the IDM, two minima can coexist. In order to guarantee the inert vacuum to be global, 
we therefore require [U 021 03] 


"^11 > ”4 


( 8 ) 



3.2. Experimental constraints 

Apart from the theoretical constraints discussed above, several experimental constraints exist 
which impose further bounds on the models parameter space. 

• We set the mass of the SM-like Higgs boson h to 

Mh = 125.1 GeV, 

in agreement with the results from the LHC experiments [33] = 125.09 ± 0.24 GeV. 

• We furthermore require the total width of the 125 GeV Higgs to obey an upper limit |5ll|52| 

Ttot ^ 22MeV. 

The validity of the above constraint relies on the assumption that no new particles enter in 
the loop-induced ggh vertex and the corresponding continuum contributions®; both of these 
assumptions are fulfilled in the IDM. Obviously, another contamination of the proposal in |53| 
would be additional new physics diagrams leading to the same final states; in our case, this 
could come from e.g. A A production and successive decays with negligible missing transverse 
energy. As discussed below, contributions from these final states are usnally relatively small, 
and we therefore consider these to be negligible. A more detailed discussion is in the line of 
future work. 

• Furthermore, we take into account strong bounds provided by the total width of the elec- 
troweak gauge bosons (cf. e.g. |55|), in a following simple form: 


Ma,h + > mw, Ma + Mh > raz, 2Mff± > mz- 


(9) 


® The latter is implemented as a standard constraint in | 50| . 

® See the original proposal in |5.1) as well as the discussion in [54]. 
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3 Constraints 


3.2 Experimental constraints 


• We also require a 2cj (i.e. 95% C.L.) agreement with electroweak precision observables, 
parametrized through the electroweak oblique parameters S, T, U [56H59] . 

The constraints from the electroweak oblique parameters S, T and U are included by eval¬ 
uating 

xItu = x'^C“^x, (10) 

with = {S — S,T — T,U — U). The observed parameters are given by [nU] 

5 = 0.05, f = 0.09, t/ = 0.01, (11) 


while the unhatted quantities denote the model predictions, which we obtained from inter¬ 
facing with the publicly available code 2HDMC [50]. The covariance matrix reads |60j ^ 


(C),y 


/ 0.0121 0.0129 -0.0071 

0.0129 0.0169 -0.0119 

-0.0071 -0.0119 0.0121 


( 12 ) 


We then require xl^u — 8.025, corresponding to a maximal 2a deviation given the three 
degrees of freedom. 


• In order to evade bounds from long-lived charged particle searches, we conservatively set 
an upper limit on the charged scalar lifetime of r < 10“^ s, to guarantee decay within the 
detector. This translates to an upper bound on the total decay width of the charged scalar 
of Ttot > 6.58 X 10-18 GeV. 


• A bound on the lower mass of Mf{± has been derived in |63) . Although this bound does lack 
a dedicated analysis within the current models’ framework, we take Mfj± > 70 GeV as a 
conservative lower limit. 


• We also require agreement with the null-searches from the LEP, Tevatron, and LHG experi¬ 
ments. We do this via the publicly available tool HiggsBounds-4.2.1 |641l66j . which includes 
results from all relevant experimental searches until summer 2014. 

• We require agreement within 2 a for the 125 GeV Higgs signal strength measurements. 
For this, we make use of the publicly available tool HiggsSignals-1.4.0 m, and require 

< 11.3139, corresponding to a 95% confidence level®. 

• Finally, some collider searches for dark matter, e.g. in supersymmetric models, could be 
reinterpreted within our model®. While doing this is beyond the scope of the current work. 


^ We here use the best linear unbiased estimator, cf. m- See also [B5]. 

* We used the signal strength parametrization of the latestresults-1.3.0-LHCinclusive data sample within 
HiggsSignals . 

® Several tools exist which allow for recasting experimental results in such frameworks, cf. e.g. |68[I69| . 
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3 Constraints 


3.2 Experimental constraints 



FIG. 1: Upper limit on dark matter-nucleon scattering cross-sections from the LUX experiment as a 
function of the dark matter mass. Taken from m- 

we use one important reinterpretation of a LEP analysis [TOj within the IDM framework |43] . 
This particularly rules out all regions where 

Ma < 100 GeV, Mh < 80 GeV, AM(A, H) > 8 GeV, (13) 


simultaneously. 

• After taking into account all the above limits we are outside of the region excluded due to 
the recent reinterpretation of the supersymmetry analysis from LHC run I [ 71 ]. 

• We apply dark matter relic density limits obtained by the Planck experiment |72|: 

= 0.1197 ± 0.0022. (14) 

In this work, we do not require the dark matter candidate of the IDM to provide the full 
relic density, but use it as an upper limit^*^ (cf. e.g. also the discussion in |20|). Being 
conservative, we require 


< 0.1241, (15) 

which corresponds to not overdosing the universe at ~ 95 % confidence level. Additionally, 
we identify regions which are within a 2 u band from the central value and thus reproduce 
the observed DM density at 95% confidence level. 

• To respect the most stringent direct detection limits from dark matter-nucleon scattering 
provided by the LUX experiment, (see Fig. 1 IZD), we use an approximation in form of an 


In such a scenario, additional dark matter candidates would be needed in order to account for the missing relic 
density; cf. e.g. m for a dedicated discussion of such scenarios within a supersymmetric setup. 
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3 Constraints 


3.2 Experimental constraints 


Constraint 
Higgs width of Mh 

no electroweak gauge boson decay into new states 

S,T,U 

decaying 

lower bound on M^± 

agreement with collider searches 

agreement with signal strength measurement 

recast LEP searches 

recast LHC searches 

agreement with dark matter relic density (upper limit) 
agreement with LUX upper bound on direct detection 


Criterium 
Th < 22MeV 
Eqn.(|^ 

< 8.025, assuming 3 d.o.f 


Used Tools 
2HDMC 
2HDMC 
2HDMC 


Ejja > 6.58 X 10-i®GeV 


2HDMC/ Madgraph 


Mh± > 70GeV 


see 


< 11.3139, assuming 5 d.o.f. 


Eqn.(13) 


not applicable after other constraints 
Vlh3 < 0.1241 (2 a deviation) 

Eqn. (16) 


HiggsBounds 

HiggsSignals 


MicrOmegas 

MicrOmegas 


TABLE I: Experimental constraints, exclusion criterium applied, and tools used for the calculation of the 
respective criterium for the experimental constraints discussed in Section 3.2 No correlation was assumed 
between the different rows listed above. 


analytic formula 
{Mh) 


cr, 


cm^ 


= 1.86934 X 10 

+ 


-46 


+ 1.1236 X 10"^^ Mn/GeV + (-1.32312 x 10“^^) (Mjy/GeV) 


1.85167 x 10"28 


(-9.35439 X IQi'^) (M^/GeV) + 1.12151 x IQi'^ (Mjy/GeV)^’ 
which, for Mh > 10 GeV, reproduces the actual exclusion limits on the 1 — 8 % level. 


(16) 


The above exclusion constraints, together with the respective constraint applied for each criterium, 
are again summarized in Tablejl} The statistical treatment of search constraints from limits provided 
by HiggsBounds and HiggsSignals have been described at great length in |641l67j . which we refer 
the reader to for more details. All above constraints are applied assuming no correlation between 
the different limits. 

Note that, for a multi-component dark matter scenario, the above upper limit depends on the 
actual DM relic density for the specific point in parameter space; therefore, we have to introduce 
a rescaling factor, leading to the (relic density dependent) limit 

oPlanck 

.(M„)<.“bM«)x ( 17 ) 

where we use Mh as a short-hand notation to describe the dependence on a specific parameter 
point The above rescaling, however, only applies in a multi-component dark matter scenario. 
We will therefore use the hard LUX limit in the remainder of the discussion, and comment on the 
effect of assuming a multi-component DM scenario in Section 


“ See also [73l[75H78]. 
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4 Scan setup 


3.3. Comparison with previous analyses 

Before the discussion of the scan setup and respective results, it is worthwhile to compare the list 
of constraints as given above with previous scan results from the literature. We restrict ourselves 
to work after the discovery of the Higgs boson, which set the SM-like Higgs mass ~ 125 GeV 
accordingly [a ESI [IZl [251 ig (see appendixfor a more detailed comparison). 

We summarize the major differences as follows 

• In this work, we apply the most updated bounds from collider and astroparticle physics. We 
consistently use limits from the Higgs signal strength by the means of a combined ht, as 
implemented in HiggsSignals, and combine these with the most updated constraints from 
astrophysical data. 

• We specifically identify regions where single constraints function as decisive, and point to 
direct correlations between theoretical or experimental constraints and limits in two dimen¬ 
sional parameter planes, whenever possible. This will help to directly determine the effects 
of improved measurements on the models parameter space, and contrast these to regions 
that are excluded per se due to theoretical constraints, and therefore will not be affected by 
possible future improvements in experimental precision. 

• We use dark matter relic density as an upper limit, which opens up regions with intermediate 
dark scalar masses < 600 GeV, which promise interesting from a collider perspective. 


• Finally, none of the studies provide detailed benchmark points for the current LHG run, 
which are highly requested by the experimental community in order to tune and optimize 
search strategies'^. Most previous studies of LHC phenomenology of this model, although 
presenting a generic idea, mainly discuss parameter points which are no longer in compliance 
with current experimental bounds (cf. discussion in Appendix . In addition to providing 
an explicit scan, the specihcation of benchmark points and planes for the present LHC 
respecting all present bounds run is a major outcome of the present study. 


4. SCAN SETUP 

As an input for our scan we chose the physical parameters; the dark scalars’ masses and 
couplings A 2 and A 345 , corresponding to self-interaction in the dark sector and interaction between 
the dark matter and SM-like Higgs boson, respectively: 

Mff, Ma, Mh±, X 2 , A 345 . (18) 


We thank the conveners of the Higgs Cross Section working group for encouraging us to select and present bench¬ 
mark points for our model. 
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4 Scan setup 


The mass of the SM-like Higgs particle was always fixed to the value Mh = 125.1 GeV. We let 
all other mass parameters float in the regime [0; 1 TeV], where we however enforce Mh to be the 
smallest of all dark scalars masses. We also impose a direct hard lower limit of Mh± > 70 GeV, 


as well as a minimal mass difference Mh± — Mh > 100 MeV (cf. discussion in Section 5.2) 


Furthermore, for Mh < 80 GeV, the combination of limits from electroweak gauge boson decays 
widths and direct LEP searches leads to a strongly constrained region in the {Mh, Ma) plane. 


which we also used an input in our scans (cf. discussion in Section 5.2). If not stated otherwise, 
A 2 G [0; 4.5] and A 345 G [—1.5; 2 ]. The range for A 345 chosen here is motivated by a set of pre-scans, 
where we did not find any allowed parameter points outside this range. If not stated otherwise, 
points were generated using flat distributions in parameter space. 

Our exclusion scan was performed in three steps. All points, whether allowed or excluded, are 
kept; all exclusion criteria for a specific parameter point are memorized. 


In the first step, we test all theoretical constraints as discussed in Section 3.1, as well as 
the total width of the 125 GeV Higgs h, the charged scalar, and EW gauge bosons. We 
furthermore include constraints from electroweak precision observables. All quantities are 
calculated using 2HDMG. Points which have passed all these constraints will be labeled ”OK 
step 1 ”. 

Only points which have passed the above constraints are then checked against limits from 
the Higgs collider searches and signal strength via HiggsBounds / HiggsSignals. Points 
which are also accepted by these cuts are in set ” OK step 2”. 

Points which passed all above bounds (unless specifically stated differently), are then fur¬ 
thermore tested against dark matter, i.e. Planck and LUX, data. This is done by processing 


through MicrOmegas (version 4.2.3) [79] and confronting with Eqns. (15) and (16) respec¬ 
tively. Points which pass these constraints are labeled ” OK step 3”. 

Points allowed by the above constrains (”OK step 3”) were then used to calculate their cross- 
sections for dark scalar pair production at the LHG. We refer to Section]^ for further details. 

In the following, whenever we quote hard limits on quantities, we want to emphasize that in 
each case we have performed extensive additional scans where the respective relation was explicitly 
violated. Eor constraints resulting from step 2 or 3, we have for each of these generated at least 
10^ such points after step 1. All statements concerning hard limits should be read in this spirit. 


The above numbers refer to general scan results reported in Section 5.4 


Einally, we want to comment on the general strategy of our scan. In our results, we have 
used 95%C.L. exclusion bounds whenever applicable, i.e. especially for bounds from direct collider 
searches as well as direct hard limits^^. We neglect correlations between these different bounds, 
and especially do not try to determine a best fit point within the parameter space^^. Therefore, 


For similar approaches, cf. |28[ I29 | for analyses with a singlet-extended Higgs sector, or [80] for a study within a 
supersymmetric model. 


Studies in this direction have e.g. been performed in [Hi¬ 
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5 Exclusions 


the density of points presented in our results section should be read having in mind we applied a 
flat scan procedure with the parameter bounds discussed above, and care must be taken to not 
interpret them in any probabilistic way. The main goal of our work is the identihcation of regions 
in the five-dimensional parameter space which are still viable after all current constraints are taken 
into account, but we do not aim at a probabilistic interpretation within these regions. 


5. EXCLUSIONS 

As stated in Sectionj^ the IDM has 5 free parameters and not all of the bounds can be translated 
into clear limits within the 2-dimensional exclusion plots. However, some of the constraints lead 
to direct limits on the parameter space, and in the following we briefly list those and the relevant 
physical reason for it. In the following, plots containing information about exclusion during step 
1 correspond to 10^ or more parameter points; all others contain 10^ points belonging to the ”OK 
Step 1” set. Furthermore, color coding usually refers to the constraints excluding certain points in 
parameter space; points which survive all constraints shown in a certain plot are presented in red. 
For more details, see the respective figure captions. 


5.1. Limits from theoretical constraints 

From the bounds tested in the first step of the scan, the requirements of perturbativity as well 
as positivity set relatively strong limits in the (A 2 , A 345 ) plane, as is visible from Fig. Here, the 
most obvious constraint stems from the perturbativity limit on the quartic coupling of the dark 
matter candidate H: Xhhhh = 3 A 2 , and therefore A 2 < ~ 4.19. Furthermore, positivity 

clearly forbids A 345 < — 1 in the whole region considered here. There is an additional lower bound 
on A 2 , stemming from positivity of the potential, which is dependent on the value of A 345 . 

None of the other theoretical constraints tested in step 1 of our scan, although putting strong 
limits, lead to clear distinctions in two-dimensional parameter planes. However, the step 2 and 
step 3 constraints, which will be discussed below, put such clear bounds in the {X 3 A 5 , Mjj±) and 
{Mh, A 345 ) planes, respectively. 

5.2. Limits from colliders 

At the second step of the scan, we consider bounds from direct scalar searches as well as the 125 
GeV Higgs coupling strengths, where we employ HiggsBounds and HiggsSignals, respectively. 
These lead to direct bounds on A 345 and A 3 together with or Mh- 

• The limits implied by HiggsBounds and HiggsSignals, being sensitive to the h decaying 
to two photons, are most obvious when considering the (A 345 , ) plane. In order to 

investigate this dependence in more detail, we allow A 345 to take values G [— 2 ; dvr], exceeding 
the allowed range for negative values of A 345 discussed above and keeping in mind that this 
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5 Exclusions 


5.2 Limits from colliders 



FIG. 2 : Limits from positivity and perturbativity constraints in the (A2, A345) plane. We here set A345 < 2 . 
The allowed points belong to the set ” OK step 1 ”. 



FIG. 3 : Step I Points which are ruled out by LHC direct searches and Higgs Signal strength, in the 

(A 


M„ 


) (left) and (A345, —m^) plane. The line shown here corresponds to Eqn. ( 19 ). Green points 

for small A3 or A345 values and large mass mass ratios are forbidden by the branching ratio h ^ H H, 

77 


with Mh < Mh/ 2 . Forbidden points in the green (lower) triangular shape are excluded from the h 
branching ratio. 


parameter space will be subjected to much stronger constraints once dark matter direct 
searches are taken into account. 


The observed limits can be explained by the fact that the virtual contribution to the h —?• 77 
decay stemming from the charged Higgs loop contains terms proportional to f {^h±) 

[a El]. A direct limit in Figure® given by 


Mh± 

Mh 


> 0.2 A 345 + 0.75, 


(19) 
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5 Exclusions 


5.2 Limits from colliders 


surviving points 
e.xcluded points 
current limit 



0 10 20 30 40 50 60 

M„lGeV] 



M„IGeVJ 


FIG. 4: Left: Total width of the 125 GeV Higgs h as a, function of the dark matter mass Mh, together 
with the current limit F/j < 22MeV. Right: A 345 as a function of the dark matter mass M^, with points 
which do (red) or do not (cyan) survive the condition on F/j; points which survive have IA 345 I < 0.1. Points 
with IA 345 I < 0.15 lead to maximal widths O (0.1 GeV). 

determines a region of parameters always allowed by the /i —>■ 77 branching ratio^®. The 
forbidden points in this figure, which apparently violate the above bound, correspond to 
parameter points where Mh < M/j /2 and are ruled out by limits on the invisible decays of 
the 125 GeV Higgs. 

• Also the total width of the 125 GeV Higgs boson h clearly sets strong limits as soon as Mh < 
Mfi/ 2 , what is visible from Figure]^ A relatively narrow stripe with IA 345 I < 0.05 — 0.1 
survives where the additional, ie. invisible, decay mode is suppressed. Note however that 
these limits for A 345 are completely superseded by bounds from direct DM searches, discussed 
below. 

• Another important constraint stems from the total decay width of the electroweak gauge 
bosons, which are implemented according to Eqn. (§. Here we need to distinguish three 
cases where a combination from LEP searches |43j and the decay widths of W and Z lead 
to quite narrow regions in parameter space: 

- Mh e [0;41 GeV]: Ma > 100GeV, 

- Mh e [41;45GeV]: Ma G [mz - Mh] Mh + SGeV] oi Ma > 100 GeV, 

— Mh G [45; 80 GeV]: Ma G [Mh] Mh + 8 GeV] or Ma > 100 GeV. 

• We take as a benchmark lifetime r < 10“^s to ensure decay inside the detector. In the 
IDM, this leads to constraints of the parameters entering the possible decay modes 

^ (A/H). 


See also Appendix for limit setting by > 77 only. 


14 







5 Exclusions 


5.2 Limits from colliders 


LUX • 
RelDens • 
OK 
OK++ 

LUX limit 



0 200 400 600 800 1000 

M„lGeV] 


FIG. 5: Results of step 3 scan on set ”OK step 2” points: allowed (red) and forbidden (dark red) regions in 
the {Mh, A 345 ) (mass in GeV) plane from the direct DM detection from LUX (dark red) and relic density 
from Planck measurements (black). Additionally points which reproduce the whole observed relic density 
were marked in gold (OK++). 


As the respective couplings are determined by the electroweak SM sector (cf. appendix Q , 
this basically puts constraints on the mass hierarchy. We found that this requirement rules 
out regions for which 


Ma > and Mu± — Mh < 100 MeV, (20) 

simultaneously hold. We have therefore enforced 

Mh± - Mh > 100 MeV 


throughout our scan. 


5.3. Limits from astrophysical measurements 

In fact the strongest limits on the IDM parameter space stem from requiring agreement with 
the dark matter measurements, included at step 3 in the scan. In general the direct dark-matter- 
nucleon scattering^® leaves a quite narrow allowed strip in the (A 345 , Mh) plane especially in the 
low dark matter mass region, cf. Figurej^ For this, rough estimates of upper LUX limits translated 
to the values of A 345 , valid for high mass region Mh ^ M^ 

A|j° 2^ ~ 0.03985 - 6.786 x lO""^ - 4.828 x 10"^ M|r, (21) 


See [82] for a detailed discussion of the calculation of the scattering cross-section within MicrOmegas. 

For lower Mh, the limit given by the fit is slightly too strong, however in this range strong constraints on low 
A 345 from too high relic density start to play an important role and this slight overconstraint has no impact on 
the results. 
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5 Exclusions 


5.3 Limits from astrophysical measurements 



FIG. 6: Step 1 points below 100 GeVafter the application of collider, i.e. HiggsBounds /HiggsSignals, 
and dark matter limits. Points below 60 GeV were sampled more frequently. 
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FIG. 7: DM density as a function of Mh (left) as well as relative difference of Ma and Mh as a function 
of Mh (right). The exceptional points all exhibit a strong degeneracy between Ma and Mh, leading to 
dominance of co-annihilation and thereby a smaller dark matter relic density. 

A345 = —A34'^ as CTDirDet ~ IA345P. ( 22 ) 

Constraints from dark matter relic density (cf. Eqn. ©), only start playing a role in either very 
high or very low DM mass regions. As before, the region where Mh < M/j/2, which allows for the 
invisible decay of the 125 GeV Higgs, is of special interest. Here LUX data put strong limits on 
the [Mh, A 345 ) plane, (cf. Figure [^. At step 3 of our scan we found quite strong constraints on 
the allowed parameter space for Mh < 54 GeV, cf. Fig.[^ basically ruling out all points for which 
Mh ^ 45 GeV. This can be traced back to the fact that the dark matter density rises strongly 
with smaller masses, rising above 100 at DM masses around 10 GeV, cf. Fig. (left). However, 
a couple of exceptional points prevail. For these, the masses Mg and Mh are degenerate, with 
differences 0(10%), cf. Fig. (right). In this case, the dominant channels are co-annihilation 
into qq final states, leading to a much smaller total dark matter density. The tension between the 
collider limits on the BR(/i —)> HH), allowing only for very low values of A 345 , cf. Fig. [^and relic 
density measurements is visible in this region. On the other hand, for larger dark matter masses 
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5 Exclusions 


5.3 Limits from astrophysical measurements 



FIG. 8 : Values for relic density after all other constraints have been taken into account (i.e. after step 2), 
as a function of A 345 (left) as well as Mh (right). The additional line signifies agreement with lower limit 
on the relic density at 95% C.L., specified in Eqn. (14). In the dark scalar sector, a mass degeneracy was 
enforced (see main body of text for details). 
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FIG. 9: Left: The relic density as a function of dark particle mass with the dominant production channel 
(see text for explanation). Right: Leading channel contribution in plane. All points are from 

set ”OK step 2”. Scan points as in Fig.[^ 


50 GeV < Mh < 600 GeV, relic density data poses no constraint if we do not additionally require 
agreement with the lower bound on the relic density. For higher masses, however, this measurement 
again becomes quite constraining, as shown in Figure In summary, we can therefore say the 
dark matter data poses a lower bound on the mass of the dark matter candidate of Mh ^ 45 GeV, 
which obviously induces the same lower limit for Mg and Mh± ■ This is one of the main results of 
our paper. 

To further exemplify the constraints stemming from dark matter relic density, in Fig. we 
plot the dependence of it as a function of A 345 (left) as well as the dark scalar mass Mh (right), 
for points which have passed the constraints in the second step of our scan. For Figs. and 
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5 Exclusions 


5.4 Results 
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FIG. 10: As Fig. 1^ (left), but now additionally applying limits on A 345 from the LUX direct search. In the 
parameter space allowed after all constraints, pair-annihilation into vector bosons dominates. 


allowed points in the low mass region. To better understand the major contributions to the dark 
matter density, in Fig. we plot the channel identified as dominant^® contribution for the dark 
matter density, as a function of the dark matter mass, for points which have survived after step 2 
of our scan. The respective cross-sections mainly depend on electroweak gauge couplings as well 
as A 345 . One can observe (cf. e.g. [31]) that pair-annihilation to electroweak gauge bosons as 
well as SM-like Higgses dominate in large regions of parameter space, while in the low mass range 
(co) annihilation into dd and bh final states dominate. Once the limit from LUX is applied on 
^ 345 ) the dominant channel is annihilation into vector boson final states, with a few parameter 
points corresponding to annihilation into hh final states, resulting from parameter points with 
A 345 < -0.5 and Mh > 600 GeV, cf. Fig. [^ 

Generally, dark matter relic density is well below the value measured by the Planck collabora¬ 
tion (cf. Eqn. 14). In fact, we found that the measured 95% GL value of the relic density can only 
be reproduced in either the low mass {Mh < 60 GeV) or high mass {Mh ^ 600GeV) region^®. 
However, to exactly match the Planck measurement within a 95% confidence level allowance 
requires an extreme fine tuning, and generally points with masses > 500 GeV are relatively 
challenging for collider searches. We will comment on a low scale scenario which can in principle 
reproduce He in Section]^ 


5.4. Results 


In this section, we summarize the results we have obtained from the scan as described in the 
previous sections, and provide scatter plots for allowed and excluded regions in several planes 


’’Dominant” here signifies the largest contributions; next leading channels, even if contributing in a similar range, 
are not shown. 

Similar results have been obtained in | 83 |. 
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FIG. 11: Mh vs A 345 (left)a,nd Mh± vs. A 345 ( right). All exclusion limits are shown, with the scan steps 
as described in the main body of the text. 


depending on two or three of the independent input parameters Mh, Ma, Mh±, A 345 . The order 
of plotting always corresponds to the legends read from top to bottom in the respective scatter 
plot. Furthermore, as we here use results for a generic flat scan with masses up to 1 TeV for Mh, 
some specific features of the low mass region are not prominent. For example, a general flat scan 
usually features dark masses Mh A 200 GeV. We omit A 2 , as this obeys quite strong limits from 
positivity and perturbativity (cf. Fig. but otherwise does not have an impact on quantities 
related to current collider and astrophysical measurements. 


As we have seen in the previous sections, in the general scan^° the most obvious constraints stem 
from LUX limits on the coupling A 345 . In Figs. im and| 12 [ left, we show the respective constraints 
and allowed points in the dark scalar mass versus A 345 planes. We have seen that the parameter 
space opens up for larger masses Mh- This naturally leads to a higher density of allowed points in 
the higher mass region^^. The triangular shape of allowed regions appears in Fig. |ll| which gets 


slightly distorted in the other dark scalar mass plots, cf. Figs. 11 (right) and 12 (left). 


A second major new result of our study is that we found that all dark scalar masses are quite 


degenerate for Mh A 200 GeV. This becomes apparent in Figs. 12 and 13 where we plot our 
results in the dark scalars’ mass planes^^. In contrast to the constraints on A 345 , we here could not 
pin down a single constraint which leads to this hierarchy, but rather found this to be the result 


of the interplay of several constraints. Fig. 14 shows the ratio of these masses, which we found to 
be < 2 for dark masses > 200 GeV^^. We also found that Mh± > Ma for all points in our scan 


The region where Mh < Mh/2 is special and is treated separately. 

Using a non-flat distibution in the scan parameters could in principle reduce this, but is beyond the scope of the 
present work. 

For Ma and Mh±, a similar result has already been presented in |84| . 

In the low-mass region, with Mh 100 GeV, ratios of Ma/Mh ~ 10 and Mh±/Mh ~ 10 are perfectly viable, 
while < 2.4. 
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FIG. 12: Ma vs A 345 (left) and Ma vs Mh± (right). All exclusion limits are shown, with the scan steps as 
described in the main body of the text. 
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FIG. 13: Mh vs Mh± (left) and Mh vs Ma (right) All exclusion limits are shown, with the scan steps as 
described in the main body of the text. 


after step 3, cf. e.g. Fig. 14 


Generically, the available parameter space closes rapidly for lower masses of Mh, mainly due 
to the strong constraints from the LUX experiment, cf. Fig. Therefore, a general scan does 
not find many allowed points in lower mass ranges, which can of course be cured by enforcing 10^ 
valid points after step 1 or introducing more specific scan boundaries. However, even in this case 
a very low fraction survives, cf. Fig. On the other hand, for relatively low masses the relic 
density value, as given in Eqn. (14), can be achieved. For the low-mass dark matter region, where 
Mh a Mh/2, the interplay of the 125 GeV Higgs width and constraints from relic density (cf. 
Fig. [7]) lead to an extreme finetuning. The lowest allowed mass for the DM candidate was found 
to be Mh ~ 45 GeV. A more dedicated discussion of this region is in the line of future work. 
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6 Benchmark points and planes for dark scalar pair-production at a 13 TeV LHC 
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FIG. 14: (left) and (right) vs A 345 All exclusion limits are shown, with the scan steps as 

described in the main body of the text. 

6. BENCHMARK POINTS AND PLANES FOR DARK SCALAR PAIR-PRODUCTION 
AT A 13 TEV LHC 

After having provided constraints on the models parameter space, an obvious question is 
whether this model can be tested at the current run of the LHC. Therefore, we here briefly comment 
on the collider phenomenology, concentrating on a 13 TeV proton-proton collider. In the following, 
we provide total production cross-sections for pair production of the dark scalars, which can already 
give a hint of regions which could eventually be accessible at a 13 TeV collider; however, the results 
presented here do not contain a dedicated phenomenological study, including e.g. simulation of 
SM background or kinematic distributions and cuts on the decay products of the dark scalars. 
Although this is obviously an important topic to investigate, we consider this an independent 
extension of our current study, which is in the line of future work, with the current results as input. 

The D{Z 2 ) symmetry forces dark particles to be-pair produced; therefore, production cross- 
sections for 

pp^ (23) 

are of interest here. We always obtain a pair of dark matter scalars H in the final states, so the 
collider signature is in the form of 

pp X + $T, (24) 

where X denotes either a single or several SM particle(s) and stands for missing transverse 
energy (MET). In principle, the process 

pp H H (25) 

could also be investigated; however, at leading order this would lead to a collider signature of 
missing transverse energy which requires special treatment, and we therefore disregard this process 
here. 
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6 Benchmark points and planes for dark scalar pair-production at a 13 TeV LHC 


Benchmark planes 

Taking all the results presented above into account, we have calculated the leading order pair- 
production cross-sections for dark sector scalars using MadGraphS |85] with the IDM UFO model 
file |25j . where we have additionally implemented the ggh effective vertex, which describes the 
one-loop induced ggh coupling in the mtop oo limit in a standard way^^. We generally find the 
following ranges for the cross-sections of the pair-produced particles 

< 0.03 pb, 

< 0.03 pb, 

< 0.015 pb, (26) 

< 0.01 pb, 

< 0.0015 pb. 

These processes are all s-channel mediated according to 

q q ^ Z ^ H A, 
qq^ ^ H+H-, 

q^ ^ (27) 

where q' denotes a different quark flavour. The production of H~ and A A are (additionally) 
mediated by 


pp H A 

pp ^ 

pp —)• 77^ A 

pp H- 

pp —)• AA 


(gg), (qq) h {h+ h ), {AA). 


(28) 


All processes in (27) are purely SM gauge coupling induced; the BSM dependence of the pro¬ 


duction cross-sections therefore purely stems from the masses of the produced particles via phase 


space. The processes in (28) depend on the couplings A 3 and A 345 , respectively^®; in the parameter 


range allowed after the scan, we found that this contribution is relatively small compared to the 


’SM-gauge bosons induced” production processes in (27). For the channel AA, on the other hand. 


the dominant contribution is from gluon-induced Higgs production. However, we found that the 


production cross-sections for direct pair-production for AA final states is relatively small, cf. (26), 


and therefore neglect this process in the following discussion. In Figures 15 and we show results 
for pair-production total cross-sections at a 13 TeV LHC for the dominant modes in the respective 
mass planes. They constitute the benchmark planes which should be investigated in the current 
LHC run. We can see that the main factor driving the cross-sections values are masses of dark 
particles and their hierarchy. 


This is basically an adoption from the implementation in |86| : TR wants to thank E. Maina for useful comments. 
See Appendix A for the relevant Feynman rules. 
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6 Benchmark points and planes for dark scalar pair-production at a 13 TeV LHC 
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FIG. 15: Benchmark planes: Production cross-sections in pb at a 13 TeV LHC for HA (left) and 
H~ (right), in the Mh, A 345 plane (upper low) and in the Mh, Ma or Mh, plane (lower row). 


For collider phenomenology, the decays of the unstable particles A need to be taken into 
account. In our scan 


BR (A Z H) = 1, 


while for most parameter points the decay 




W^H 


(29) 


dominates; the second decay mode is given by —)• A. For the process (29), we plot the 


respective branching ratio in the {Mh, Mh±) plane in Fig. 17 All branching ratios have been 
obtained using 2HDMC. The electroweak gauge bosons decay as in the SM. As before, the couplings 
involved all stem from the electroweak gauge sector. 

Before proceeding to the discussion of benchmark points, we want to briefly comment on the 
results reported here. As discussed above, we do not perform a dedicated phenomenological study, 
which is clearly beyond the scope of the present work; especially for SM backgrounds, dedicated 
searches at 13 TeV, including the latest higher order predictions for such processes, would be 
needed, which requires an extensive computational setup clearly beyond the scope of the present 
study. However, even from the numbers quoted above one can draw conclusions as follows 
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6 Benchmark points and planes for dark scalar pair-production at a 13 TeV LHC 
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FIG. 16: Benchmark planes: Production cross-sections in pb at a 13 TeV LHC for (left) and 

A (right), in the Mh, A 345 plane (upper low) and in the AIh plane (lower row). 
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FIG. 17: Branching ratio of charged scalar decay into W boson and dark matter (H). The only other 

channel is —>■ AW^. 

• According to the current design strategy at LHC Run H (cf. e.g. [H7]), at 13/14 TeV the 
LHC aims at a total integrated luminosity of ~ 100 fb~^ before the second long shutdown 
(LS2, in 2019) and a total integrated luminosity of about 300 fb~^ before the third LS (in 
2023/4). 
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6 Benchmark points and planes for dark scalar pair-production at a 13 TeV LHC 


• These numbers can be easily translated to number of events which could be produced for 
these two time benchmarks: for about 100 total events produced at LS2 (LS3), a minimal 
total production cross-section of 1 (0.3) fb would be needed. 

• Considering the numbers presented above, this means that a priori all channels discussed 
here can and should be investigated by the experiments at LHC run H. However, in practise, 
for LS2 only scenarios for dark masses < 400 GeV might provide large enough cross-sections 
to render feasible predictions. 

Although a more dedicated study might provide more information, we consider the above already 
an important input for the explicit design of experimental search strategies at LHC Run H. 


Benchmark points 

In addition to the benchmark planes presented above, we also provide some benchmark points 
for our model which should be considered for investigation at the present LHC run^®; other viable 
benchmark points for masses < 1 TeV have e.g. been presented in [laissiisi]. For this, we fix the 
dark scalars masses. The couplings A 2 , A 345 , when kept within the allowed ranges discussed above, 
do not have a significant influence on the production and decay processes discussed here^^. For 
these, we therefore give the allowed ranges which should be kept for consistency reasons. Cross- 
sections have been produced as described above. In the list of benchmarks given below, we provide 
respective production cross-sections for dark scalar pair-production at a 13 TeV LHC. 

• Benchmark I: low scalar mass 

Mh = 57.5 GeV, Ma = 113.0 GeV, = 123, GeV, A 2 G [0;4.2], IA 345 I G [-0.002; 0.015] 

HA : 0.371(4)pb, H+ R- : 0.097(l)pb , : 0.3071(4)pb, A : 0.1267(l)pb 

A decays 100 % to Z H, and —)• VF+ H with a BR > 0.99. 

• Benchmark II: low scalar mass 

Mh = 85.5GeV, Ma = 111.0GeV,= 140, GeV, A 2 G [0;4.2], IA 345 I G [0.;0.015] 

HA : 0.226(2)pb,i7+iL- : 0.0605(9)pb , H^H : 0.1439(2)pb, H^ A : 0.1008(l)pb 

A decays 100 % to Z H, 77+ —)• 1T+ 77(A) with a BR ~ 0.96(0.04). 

See e.g. |88H90| for benchmark points for general 2HDMs, as well as Appendix [^for previous benchmark suggestions 
in the IDM which are excluded by now. 

In this respect, the benchmarks above are also safe from constraints when vacuum stability and perturbativity are 
promoted to higher scales, as e.g. in |25 |. 
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6 Benchmark points and planes for dark scalar pair-production at a 13 TeV LHC 


• Benchmark III: intermediate scalar mass 

Mh = 128.0GeV, Ma = 134.0GeV,= 176.0, GeV, A 2 G [0;4.2], IA 345 I G [0.;0.05] 

HA: 0.0765(7)pb, H+ R- : 0.0259(3)pb, H^H : 0.04985(5)pb, A : 0.04653(5)pb; 

A decays 100 % io Z H, 77+ —^ H{A) with a BR ~ 0.66(0.34) 

• Benchmark IV: high scalar mass, mass degeneracy 

Mh = 363.0 GeV, Ma = 374.0 GeV, = 374.0 GeV, A 2 G [0;4.2], IA 345 I G [0.;0.25] 

H,A: 0.00122(l)pb, H+R- : 0.00124(l)pb, H^H : 0.001617(2)pb, H^ A : 0.001518(2)pb; 
A decays 100 % to Z H, and 77+ 100 % to W^H 

• Benchmark V: high scalar mass, no mass degeneracy 

M//= 311.0 GeV, Ma = 415.0 GeV, = 447.0 GeV, A 2 G [0;4.2], IA 345 I G [0.;0.19] 

77, A : 0.00129(l)pb, 77+ 77" : 0.000553(7)pb, 77+ 77 : 0.001402(2)pb, 77+A 

0.0008185(8)pb; 

A decays 100 % to Z 77, 77+ —)• 1V+ 77 with a BR > 0.99 


While benchmarks I and II are exceptional points in a sense that the allowed parameter space is 
extremely constrained in the low mass region, benchmarks III to V are more typical, as these parts 
of the parameter space are more highly populated. Furthermore, for scenario IV the production 
cross-sections for HA and 77+ 77“ have similar order of magnitude. Of all benchmarks above, only 
the first one is able to produce the correct 95% CL relic density (Eqn. 14). 

IDM searches at the LHC call for a detailed investigation of the best experimental strategy, as 
final states with IDM particles can have different topologies. For example, in the two lepton + 
MET channel, the emitted leptons can have significantly different kinematic properties, namely: 


• HA —)• HHZ —)• 7777/+Z with two boosted leptons, with an invariant mass around mz for 
on-shell decays 

• 77+77" —7- HHW~^W~ —> HH9vHI~ with two leptons which can be well separated. 


Studies of possible search strategies for the IDM were discussed in e.g. [ZmSlEZlEg. 

Additionally, we also show the total width of the unstable dark scalars in Eig. 18 For all points, 
F/M ~ 0.05, with F/M ~ 0.01 being a more typical value. We therefore expect that use of the 
Narrow Width Approximation is well justified in the parameter space explored in our study. 
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Full matrix elements including interference can of course be easily used in the Monte Carlo framework described 
above. 
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7 Multi-component dark matter scenario in IDM 



FIG. 18: Total widths of unstable dark particles, in GeV: A (left) and (right) in plane of their and 
dark matter masses. For most points, T/M < 0.01, with the maximal ratio F/M ^ 0.05. 
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FIG. 19: Step 1 points for Dark Matter masses between 60 GeV and 100 GeV after the application of 
collider, i.e. HiggsBounds /HiggsSignals, and dark matter limits. Points below 60 GeV were sampled 
more frequently. 

7. MULTI-COMPONENT DARK MATTER SCENARIO IN IDM 

Before concluding, we want to briefly comment on changes in our analysis when we consider 
multi-component dark matter. In this case, the very strict limits from LUX in the (Mh, A 345 ) are 
significantly altered, and the allowed parameter space is greatly enhanced. This concerns especially 
upper values for the A 345 coupling. Below we plot several of our previous figures, now using the 
maximal allowed value for direct detection cross-sections via rescaling (cf. Eq. We can 

summarize our findings as follows 

• for Mh < 60 GeV, the picture basically does not change. Here, the strongest constraints 
on the {Mh, A 345 ) plane stem from the Higgs signal strength measurements, as well as relic 
density, the impact of LUX limits is minimal. However, for masses Mh £ [50; 60]GeV the 
largest values for IA 345 I are slightly enhanced, from 0.02 to 0.03 as maximal values. Note 
also that this is the region where in addition ~ 1 . 
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FIG. 20: (Left): As Fig. 11 (right); (Right): as Fig. 12 (right), with rescaled LUX limits. 


for masses G [60; 100]GeV, the picture dramatically changes. Now, A 345 can become as 
large as ~ 8 for masses in the co-annihilation regime, where the dark matter relic density is 


relatively low, cf. fig 19 Furthermore, for masses Mh ^ 80 GeV again the parameter space 
opens up significantly. 

for Mh > 100 GeV, the parameter space in the (Mh, A 345 ) plane also largely opens up, 
allowing A 345 to reach values ~ dvr, cf. Fig. The upper limit on this value are now 
basically set by S, T, U, perturbativity, and unitarity constraints. In particular, constraints 
from direct detection do no longer play a prominent role in the determination of allowed 
regions in the (Mh, A 345 ) parameter plane. 

Similarly, the relation Mh± > no longer holds; however, we still observe the same 
degeneracies between the dark scalar masses as before. 

Regarding total widths of the dark scalars, the general behaviour also does not significantly 


change, and total widths are similar to the case without rescaling (see Fig. 18). The ratio 


T /M can now reach up to 10 % for the second neutral scalar A. As before, the charged scalar 
dominantly decays into W~^ H. For the decay of A, the additional channel also 


opens up, cf. Fig. 21 


A primary effect of opening up the additional parameter space for A 345 is the fact that now the 
production channel A A also renders cross-section values which are of similar order of magnitude 
as the other channels, cf. Fig. 21 Therefore, a discovery in this channel could additionally 


give information about the nature of the dark matter scenario, which again highlights the strong 
connection between collider and astroparticle physics in the IDM. All other production cross- 
sections are mainly gauge boson induced, and the production cross-sections remain at a similar 
order of magnitude. 

Furthermore, if we now investigate the dominant decay channels responsible for relic density, 
the relaxed constraints now also allow regions where channels other than the diboson annihilation 
dominate, and nearly all points from Fig. (left) are now allowed by the LUX limit. 
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8 Conclusions 
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FIG. 21: (Left): Branching ratio of A to Z H; the other decay mode is A —)■ (Right): Production 

cross-sections in pb at a 13 TeV LHC for in the Ma, A 345 , when using rescaled LUX bounds. 


8. CONCLUSIONS 

In this work, we have revised the current theoretical as well as experimental constraints on a 
two Higgs doublet model with an exact Z 2 symmetry (Inert Doublet Model), which extends the 
SM Higgs sector in a straightforward way and in addition comes with a dark matter candidate. 
We have included all experimental bounds from direct searches, the 125 GeV Higgs coupling 
strength, as well as limits from recent astrophysical measurements such as limits from direct DM 
searches as well as dark matter relic density, the latter in the form of an upper limit. We have 
closely investigated the interplay of these bounds for dark scalar masses which are in the discovery 
regime of current and future collider experiments. We found that an interplay of astroparticle 
and collider constraints rules out regions with Mjy < 45 GeV, thereby setting a minimal mass 
scale on this model. Furthermore, the parameter space is highly constrained by direct DM search 
experiments, which supply stringent limits especially in the low mass region. 

Apart from clear cuts on the additional couplings, we found that the combination of all 
constraints leads to a mass hierarchy in the dark scalar sector. The combination of different 
constrains ensures Mff± > Mg > Mh- In addition, in the flat scan more degenerate masses of 
dark scalars are preferred for Mu > 200 GeV. Furthermore, as we do not require the IDM to 
provide the full relic density, intermediate masses of the dark matter candidate (~ 100 — 600 GeV) 
become available, which open up appealing parameter regions for the current LHG run. Our results 
therefore correspond to the most up-to-date investigation of the models parameter space, thereby 
providing important information about viable parameter regions for collider phenomenology. 
Furthermore, we also briefly commented on the case of multi-component dark matter. In this 
case, the limits from direct detection are relaxed and we find that especially for Mh > Mh/2, the 
parameter space opens up. This in turn leads to cross-sections for A A pair-production which are 
of similar magnitude as production, and could therefore in principle also be investigated 

at the current LHC run. This again emphasizes the complementarity of collider searches and 
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astrophysical measurements within the IDM. 

In addition, we have provided benchmark planes and cross-sections for pair-production of the 
dark scalars for the dominant production modes within our model in terms of the relevant new 
physics parameters, which should be used as input for the experimental searches. As both pro¬ 
duction and decay modes for the new particles are mainly determined by kinematics, the above 
degeneracies directly induce relatively small production cross-sections as the masses increase, lead¬ 
ing to (T ~ I fb for masses > 300 GeV. Nevertheless, this model constitutes an attractive extension 
of the SM, and provides an excellent example of complementarity between astroparticle and collider 
experiments. We strongly encourage the LHC experiments to investigate the benchmarks proposed 
here in the current LHC run. 
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Appendix A: IDM Feynman rules and other relations 


The parameters r?T, 22 , A 3 , A 4 , A 5 can be re-expressed in terms of our input parameters; 


^22 ~ A345n^ —2117^, 


A3 — A345-2 ~ ^ H ±) ) -^4 — 




) A5 — 


Mjj-Ml 


(Al) 

(A2) 


Eor completeness, we list the relevant Eeynman rules including dark scalars (omitting goldstone 
modes, as we are working in the unitary gauge and at tree level), see Tables [H| and HI . Note that, 
since the second doublet does not participate in electroweak symmetry breaking and fermion mass 
generation, the couplings of the SM-like Higgs h to electroweak gauge bosons as well as fermions 
are given by their SM values, see e.g. [93], with the convention = ie^v/2sw‘^9tiv • 
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vertex 

coupling 

hHH 

A345 ^ 

hAA 

A345 V 

hhh 

3 Ai V 

hH+ H~ 

X3V 

hhhh 

3 Ai 

H+ H+ H- H- 

2 A2 

HHAA 

A2 

HHHH 

3 A2 

AAAA 

3 A2 

H+H- AA 

A2 

H+H-HH 

A2 

hhH+H- 

^3 

hhHH 

A345 

hhAA 

A345 


TABLE II: Non SM-like couplings in the scalar sector of the IDM; a global factor of —i is omitted. 


vertex 

coupling 

H- 77+7 

i e 

77- H+ Z 

„■ g cos {2ew) 

2 cos 6w 

77 77+ 1T+ 

Tff 

A 77+ 1T+ 

9 

2 

HAZ 

9 

2 cos 9 w 


TABLE III: Non-SM couplings to gauge bosons in the IDM; all vertices are proportional to a kinematic 
factor (p2 — Pi)/i for a vertex of the form Ai(pi) X 2 {p 2 ) ^3(^3). 


Appendix B: as a standalone variable 

As an alternative to the generic fit using HiggsSignals, and also as a consistency check, we can 
regard the branching ratio of the SM-like Higgs into two photons as a single constraining variable. 
We then take as a reference value |94j 

= (2.28 ± 0.11) X lO-l 

Furthermore, we use m 

= 1.17 ± 0.28, = 1.12 ± 0.24 

which yields a naively combined value of 

= 1.15 ± 0.18. 
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1 Previous scans 


Again combining these, we have 

= (2.62 ± 0.43) x 10"^ 

which at 2 it yields 

e [1.76; 3.84] x 10“^. 

Obviously all errors combinations above neglect both correlations as well as non-gaussianity of the 
respective distributions. 


Appendix C: Previous scans and benchmarks 
1. Previous scans 

In this subsection, we give a brief overview on similar studies which have been performed after 
the discovery of a 125 GeV Higgs, and point to differences with respect to the constraints and 
results achieved in this work. We specifically discuss [nmni [13 nails]. 

• Obviously, all above studies respect the theoretical bounds on the potential, i.e. positivity, 
as well as constraints from perturbative unitarity. Most of them also include the constraint 
in Eqn. Q. As we did not find points in parameter space where this is the most important 
constraint, we can state that theoretical bounds where treated on equal footing^®. 

• Similarly, decays from electroweak gauge bosons, electroweak precision constraints in terms 
of oblique parameter (although with slightly different central values), a lower bound of the 
charged Higgs mass ~ 70 GeV, as well as recasted LEP results were also included in all 
above studies. 

• Refs. |14[ [T71142] particularly discuss the interplay of constraints from the 125 GeV Higgs 
signal strength measurement, mostly in terms of possible modifications to the decay rate 
to diphoton final states, with constraints from dark matter detection. In fact, these will 
also prove important especially in the low mass region of our scan, and have therefore been 
correctly identified in the above work as the most important constraints for certain parts of 
parameter space. 

• More specifically, in ref. [13112] the effect of including the by that time newly measured 
Higgs mass is discussed, and predictions for the SM-like Higgs decay into diphoton or 
7 Z final states are considered, leading to additional constraints. The authors present their 
results in terms of two-dimensional planes, and determine allowed/ forbidden regions in 


In [^, scale dependent limits where also tested at larger scales. However, if we consider the IDM as an effective 
model at the electroweak scale, the limits coincide. 
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1 Previous scans 


parameter space from the constraints included. Obviously, a more detailed analysis using 
Higgs coupling strength measurements and LUX data was not possible at that time, so these 
have not been included in the above study. 

• In ref. HZ!, the above studies are combined with constraints from XENON [96| direct de¬ 
tection limits. In this work, the authors concentrate on low or middle regions for the dark 
matter mass, and mainly consider the case where Mh < 100 GeV. Furthermore, they con¬ 
centrate on identifying regions where the IDM constitutes the full dark matter content of the 
universe. Our study clearly differs insofar as we allow for lower relic density values, which in 
turn opens up the parameter space for higher dark matter masses. In addition, we provide 
results in terms of a full random scan, which is also not done in this work. 


• Ref. [25] is closely related to our work, where the authors in addition consider vacuum stabil¬ 
ity and perturbativity at arbitrary scales using renormalization group equations (RGBs). As 
expected, requiring validity at higher scales constrains the allowed parameter space consid¬ 
erably. However, we find that after all constraints have been taken into account, parameters 
which are important for collider phenomenology are in similar ranges^®. The direct detec¬ 
tion limits applied here follow from XENON [96]. Furthermore, out of the benchmark points 
presented here only one survives the updated constraints presented in our study. 


Ref. [16] seems to be closely related to our study. In this work, the authors basically include 
all above constraints apart from perturbativity of the Higgs self couplings, the total decay 
width of the 125 GeV Higgs and the charged scalar; the inclusion of direct search bounds 
was persecuted on a slightly different footing. Furthermore, limits from signal strength 
measurements were not done on the level of a global fit as implemented in HiggsSignals. 
However, here again the authors require exact agreement with the dark matter relic density, 
while we include the latter requiring an upper bound only which immediately opens up larger 
parts of parameter space. The main focus of the authors is a presentation of the successive 
combination of constraints on profile likelihood fits, which they discuss in detail in their work. 
In our study, this is performed in a slightly different way, allowing the possibility to present 
results in terms of several individual two dimensional parameter planes directly showing the 
impact of separate constraints (cf. section 5.4). Whenever applicable, the regions they find 
consistent with all their bounds are indeed similar to the results obtained here, which aligns 
with our results. Furthermore, the authors do not provide predictions for production cross 
sections at the current LHG run, which we consider imminent. 


The only difference is that we allow for slightly larger negative values of A345 for dark masses > 600 GeV, and a 
larger range of A2. As discussed below, these parameters are of no relevance concerning LHC phenomenology. 
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2. Previous benchmarks in the literature and exclusions 

In this section, we present benchmarks which have been discussed in the literature for the IDM 
in the context of LHC phenomenology, with Higgs masses roughly ~ 125 GeV (allowing for a 
range of ±5 GeV or similar), and the reasons why these are by now excluded, following the step 
of constraints as discussed in the main body of the manuscript. We do not claim this to be an 
exhaustive list of all benchmarks ever presented. We furthermore restrict ourselves to scalar masses 
< 1 TeV, in order to maintain cross sections which might be probed at current colliders. 

• Benchmarks ”LH2, LH3, LH4, LH5” from [22] 

LH2: excluded by S,T,U, and total width of h 
all others: excluded from ATLAS search HZl 

• additional benchmarks ”LH6, LH8” from |92| 

LH6: excluded by total width of h 

LH8: excluded from ATLAS search [97] 

• IDM benchmarks suggested in |98j 
BM 1,3 : from STU 

all others: from direct detection 

• benchmarks from m 

Table H: first benchmark excluded from combination of LEP recast and decay width of 
electroweak gauge bosons; second benchmark excluded from dark matter relic density^^, 
benchmark 3 from signal strength (too large discrepancy in Mh), benchmark 4 from pertur- 
bativity 

Table HI: first two benchmarks from perturbativity, last from STU 

• benchmarks from [2] 

all ruled out from combination of LEP recast and electroweak gauge boson decay widths 

• benchmarks from [25] 

BM I: excluded from direct detection 
BM HI: excluded from STU 
BM H allowed 

• benchmark IV from |36] 

excluded by combination of LEP recast and electroweak gauge boson decays 

• benchmarks from table VH in m 
Benchmarks 1,6: excluded from direct detection 


Using micromegas 4.2.3, with exactly the same input values 
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2 Previous benchmarks in the literature and exclusions 



a■AH[p^ 

[pb] 

BR {H+ W+ H) 

2 

0.1275 (1) 

0.03949(4) 

1.0 

3 

0.9208 (5) 

0.1063 (2) 

0.68 

4 

0.2673 (3) 

0.04908 (1) 

1.0 

H 

2.1 X 10"'* 

8.6 X 10"® 

0.81 

1 

1.7 X 10"'* 

1.8 X 10-4 

0.91 

2 

7.2 X 10"® 

7.7 X 10-® 

1.0 


TABLE IV: production cross sections for valid benchmarks from [15] (first three entries), |25| (middle entry, 
with roles of A and H interchanged), [HI] (last two entries) with an adjusted mass = 125.1 GeV 


Benchmark 2: excluded from STU 
all others allowed 

For completeness, we provide production cross sections and dominant decay modes for scenarios 
which have been discussed above in table lYl 
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